Maxwell's Multipole Vectors and the CMB 
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15 Farmer Street, Canton NY 13617 US^ 

The recently re-discovered multipole vector approach to understanding the harmonic decomposition 
of the cosmic microwave background traces its roots to Maxwell's Treatise on Electricity and Mag- 
netism. Taking Maxwell's directional derivative approach as a starting point, the present article 
develops a fast algorithm for computing multipole vectors, with an exposition that is both simpler 
and better motivated than in the author's previous work. 

Tests show the resulting algorithm, coded up as a Mathematica notebook, to be both fast and 
robust. This article serves to announce the free availability of the code. 

The algorithm is then applied to the much discussed anomalies in the low-^ CMB modes, with 
sobering results. Simulations show the quadrupole-octopole alignment to be unusual at the 98.7% 
level, corroborating earlier estimates of a l-in-60 alignment while showing recent reports of 1-in-lOOO 
(using non-unit-length normal vectors) to have unintentionally relied on the near orthogonality of 
the quadrupole vectors in one particular data set. The alignment of the quadrupole and octopole 
vectors with the ecliptic plane is confirmed at better than the 2a level. 

PACS numbers: 98.80.-k 



I. INTRODUCTION 

The Wilkinson Microwave Anisotropy Probe (WMAP) 
provides an unprecedented view of tfie Cosmic Microwave 
Background (CMB) sky. WMAP's first-year observa- 
tions m, while generally matching researchers' expecta- 
tions in the high-^ portion of the CMB's spherical har- 
monic decomposition, yield suspicious results in the low-^ 
portion of the spectrum. Almost immediately after the 
first-year data release, Tegmark and colleagues noticed 
that the reported quadrupole {£=2) and octopole (^=3) 
align 0, and Eriksen et al. found statistical inconsis- 
tencies between the northern and southern ecliptic hemi- 
spheres . Eriksen and colleagues later found that while 
the i = A component is generic, the i = 5 component is 
unusually spherical at the 3<t level while the £ = 6 com- 
ponent is unusually planar at the 2a level Q. 

To make sense of these anomalies in a systematic way, 
Copi et al. introduced "a novel representation of cosmic 
microwave anisotropy maps, where each multipole order 
£ is represented by £ unit vectors pointing in directions 
on the sky and an overall magnitude" jSj] . The beauty of 
this scheme is that the CMB alone determines the multi- 
pole vectors, without reference to the coordinate system. 
This contrasts sharply with the more common aim repre- 
sentation of spherical harmonics, which depends strongly 
on the coordinate system. The coordinate-free, geometri- 
cal nature of the multipole vectors makes it trivially easy 
to test for alignments between modes (for example the 
quadrupole and octopole align with each other) and with 
external reference points (the quadrupole and octopole 
align with the CMB dipole) 0. 

Inspired by these results, Katz and Weeks went on to 
recast the traceless tensor methods of Refs. 0, into 



the language of homogeneous harmonic polynomials, and 
within that context to prove the existence and uniqueness 
of the multipole vectors and to devise an algorithm for 
computing them |2j . 

Even though Copi et al. fancied themselves the dis- 
covers of multipole vectors and Katz- Weeks fancied 
themselves to be breaking new ground in the polynomial 
approach, the truth is much older. The true discover of 
multipole vectors was James Clerk Maxwell in his 1873 
Treatise on Electricity and Magnetism ! Maxwell's 
starting point was, of course, electromagnetism. He 
began with the ^ potential surrounding a point charge (a 
monopole) and asked what happens when two monopoles 
of opposite sign get pushed together. Elementary reason- 
ing shows the potential of the resulting dipole to be given 
(up to rescaling) by the directional derivative Vuj-^, 
where Ui is the direction along which the two monopoles 
approach one another. Similarly, pushing together two 
such dipoles with opposite signs gives (up to rescaling) 
a quadrupole with potential VuaVui^, where U2 is the 
direction along which the dipoles approach, and so 
on. Maxwell recognized these potentials as solutions of 
Laplace's equation, i.e. spherical harmonics. Because 



his expression AV^ 



contains 2£ + 1 degrees 
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of freedom, and because a simple dimension counting 
argument indicates that a generic harmonic polynomial 
of degree £ also has 2£ + 1 degrees of freedom. Maxwell 
concluded - with less than perfect rigor but nevertheless 
correctly - that his method yields all spherical harmonics. 

Section III Al establishes, by purely elementary means, 
the equivalence between Maxwell's multipole vectors and 
the vectors used in more recent work ^^Jj . Section FlI Bl 
sheds fresh light on the fast algorithm for computing mul- 
tipole vectors, providing a simple, natural motivation for 
a construction that felt unnatural and ad hoc in Ref. • 
Section lTTTI announces the release of a Mathematica pack- 
age implementing the fast algorithm, and analyzes its ef- 
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ficiency and stability. Finally, Section lTVl applies the fast 
algorithm to the first-year WMAP data, finding both the 
quadrupole-octopole alignment and the alignments with 
the ecliptic plane to be significant only at the 99% level, 
not at the higher levels claimed in Refs. ^J}- 



II. THEORY 
A. Proof of equivalence 

While Maxwell's directional derivatives provide a wel- 
come breath of fresh air for researchers toiling with more 
difficult approaches^namely traceless tensors and fac- 
tored polynomials one must nevertheless prove that 
Maxwell's multipole vectors are indeed the same vectors 
that appear in those other approaches. Mark Dennis has 
successfully applied Fourier methods 0, Appendix A] 
to prove the equivalence of the polynomial and Maxwell 
interpretations. Here we obtain the same result by more 
elementary means, keeping to the spirit of Maxwell's orig- 
inal work PI Chap. IX]. 

Maxwell expresses a spherical harmonic as 



to 



fi{x, y, z) = Vuj • • • VuaVui -, 

r 



(1) 



where r = + y"^ ^ . Observe the simple pattern 
that emerges as we apply the directional derivatives one 
at a time, 



(-1)(U1T) 

(3-l)(ui-r)(u2-r)+r2(-Ui-U2) 
(-5-3-l)(ui-r)(u2-r)(u3.r)+r^(...) 

(2) 

where boldface r = {x,y,z) while plain r = ^/r ■ r as 
before. The ellipsis (...) marks a polynomial whose 
form docs not interest us. 

Let Pi denote the polynomial in the numerator of each 
ff. The action of Vu^, may be written explicitly as 
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~(2£- l)(u^ ■ r)P^_i + r2(u, ■ VF,_i) 



(3) 



It is now obvious at a glance that formula takes 

(-1)^-1 (2^ - 3)!! (ui • r) • • • (u,_i • r) + r^{. . . ) 



fe-i = 



„2£-l 



fi 



i-iy (2£~1)!! (ui-r).-.(u,.r)+r2(...) 



(5) 



thus establishing by induction the validity of the latter 
for all i. 

We would now like to understand more deeply the 
relationship between the rational function fi and the 
polynomial Pi appearing in its numerator. Following 
the method of electrical inversion, which Maxwell [3, 
art. 162] credits to Thomson (Lord Kelvin) and Tait [Tj 
but Axler [l^] traces back to Kelvin's work two decades 
earlier p^, define the Kelvin transform f of a function 
/ to be 



f{x,y,z) = -/(4- 



y_ 



(6) 



The Kelvin transform in effect reflects a function across 
the unit sphere and then adjusts its amplitude to keep 
it harmonic (to be proved below). 

Lemma 1. The functions / and / agree on the unit 
sphere r = 1. 

Lemma 2. The Kelvin transformation ^ is an involution, 
that is. 



f{x,y,z) 



y 



4) 

-^f ^!{x,y,z). (7) 

r _9 — T — T — T 



Lemma 3. If / is harmonic, then so is /. Maxwell calls 
this the theorem of electrical inversion P, art. 129]. It is 
easily proved by a direct calculation 



V' f 5' i)) = i(vv)( J, ^, ^) = 0. (8) 



r r" r" 



The converse, that / harmonic implies / harmonic, is 
automatically true by Lemma 2. 

Each multipole function fi = ^^t+i has a polynomial 
numerator P^ of homogeneous degree ^, so its Kelvin 
transform is easy to compute. 



1 P'- 

~ 1 —97- 



(9) 



(4) 



In other words, the Kelvin transformation of ft extracts 
the numerator Pi. Lemma 3 now implies that Pi is a 
harmonic polynomial in its own right. 
Equation ||SJ) gives Pi explicitly as 

Pi = i-lY {2£ - 1)!! (ui • r) • • • (u, ■r)+r'{...). (10) 

Happily this is exactly the "factored form" whose 
existence and uniqueness were established in Ref. and 
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elsewhere. Thus we have proven 
Proposition 4- Each harmonic function 

Vu, •••Vu.Vui ^ (11) 

defined by Maxwell agrees on the unit sphere with the 
"factored form" 

A(uit)---(u, •r)+r2Q (12) 

of Ref. In particular, a given function employs 

the same set of multipole vectors {ui, . . . ,Ui} in both 
representations . 

The only remaining question is whether the mapping 
from Maxwell's form Hll(l to the factored form H12() is 
surjective. For each choice of {ui, . . . , u^} the mapping 
hits exactly one factored polynomial 

A(uiT)---(uf •r)+r2Q (13) 

while infinitely many other polynomials 

A(ui • r) • • • (u, • r) + r^Q' (14) 

with Q' ^ Q, will not get hit. So strictly speaking the 
mapping is far from surjective. The reason is that the fac- 
tored form applies to all homogeneous polynomials, while 
Maxwell's method generates only homogeneous harmonic 
polynomials. If we restrict the factored form from all ho- 
mogeneous polynomials to only the harmonic ones, then 
the mapping becomes surjective for the following rea- 
son. The polynomials (|13|l and (|14|l cannot both be har- 
monic, because then their difference r^{Q~Q') would be 
harmonic as well, contradicting the fact that no nonzero 
multiple of is harmonic (Endnote 17]). Equation H10() 
gives a value of Q for which (|13|l is indeed harmonic, thus 
proving that the mapping from Maxwell's form l|ll|) to 
the factored form H12|l has as its image exactly the set of 
harmonic polynomials. 

B. The algorithm 

The most efficient algorithm for computing a spheri- 
cal harmonic's multipole vectors expresses a given spher- 
ical harmonic as a homogeneous harmonic polynomial 
P{x, y, z) and then looks at the intersection of the graphs 
P(a;, y,z) = Q and a;^ -l-y^-l-z^ = in the complex projec- 
tive plane. The Katz- Weeks paper simply pulls this 
approach out of thin air, because at the time my coau- 
thor and I did not understand the origins of our work. 
It turns out that the key ideas in "our" algorithm first 
appeared in 1876 in J.J. Sylvester's "Note on Spherical 
Harmonics" [3. While Sylvester makes his harsh feel- 
ings for Maxwell abundantly clear, he leaves his algo- 
rithm and the motivation behind it decidedly unclear. 
Inspired by Sylvester's paper, the following paragraphs 



lay out the algorithm in elementary terms, speculating 
on what Sylvester's reasoning might have been. 

Sylvester's starting point was Maxwell's directional 
derivative formulation of a spherical harmonic, which 
reduces to a polynomial (|10|1 . which for present purposes 
may be written as 

P{x,y,z) = A[|(u,-(x,y,z)) + (2;2+J/2 + 22).g^ (^5) 

1=1 

Let us take stock of the computational task that lies be- 
fore us. We are given the harmonic polynomial P{x, y, z) 
of homogeneous degree £, and we are asked to compute 
the linear factors Ui-{x, y, z), but with no prior knowledge 
of the remainder term Q. 

If we could find a common solution to P(x,y,z) = 
and -f + = 0, then such a solution would also 
satisfy at least one of the linear equations • (a;, y, z) = 0. 
Of course x^-l-y^ + z^ = has no nontrivial solutions over 
the real numbers, which motivates us to consider complex 
solutions instead. Postponing the computational details, 
assume for the moment that we have found a common 
solution {x, y, z) to P{x, y,z) = and x'^ + y"^ + z"^ = 0, 
which is therefore a solution to • {x, y, z) = for some 
i. Moreover, because is real, the real and imaginary 
parts of {x,y,z) individually satisfy the condition 

Ui ■ (Re a;. Key, Rcz) = 

Ui ■ (Imi, Imy, Imz) ~ 0. (16) 

Furthermore, the real vectors (Rex, Key, Rez) and 
(Imx, Imy, Imz) must be linearly independent, because 
otherwise {x,y,z) could not be a nontrivial solution of 
x"^ + y'^ + z"^ = 0. Therefore their cross product, when 
normalized to unit length, gives the desired multipole 
vector 

Ui ^ (Rex, Rey, Rez) x (Imx, Imy, Imz). (17) 

If we can find enough distinct solutions {x,y, z), we will 
get all the multipole vectors u^. 

So how may one find the common solutions {x,y,z) 
to P{x,y,z) = and x'^ + y'^ + z'^ = ? Sylvester 
applied his own resolving equation method, and Math- 
ematica might be employing a similar algorithm in its 
proprietary NSolve function. Unfortunately NSolve is 
slow. If we give Mathematica a little help we can speed 
the computation up enormously. Following Ref. "T], we 
parameterize the graph of -I- y^ -I- = 0, which topo- 
logically is a Riemann sphere CP^ w S^ sitting inside 
CP^, via the parameterization a {x{a),y{a), z{a)) = 
(a^ - 1, 2a, i {a^ -I- 1) ), for a G C U {00}. Evaluating 
P on the image of that parameterization yields a single 
equation P{a'^ — 1, 2a, i (a^ + 1) ) = of degree 2£ in 
the single variable a, which Mathematica solves quickly 
and accurately for the 2i values of a (counting multiplic- 
ities). Applying the parameterization gives 2£ solutions 
{xi, {ji, Zi) = ( af — 1, 2ai, i {af + 1) ). Because the orig- 
inal equations are real, whenever a point {xi, jji, Zi) sat- 
isfies them, the complex conjugate [xi, iji, Zi) will satisfy 
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them as well. Th us the 21 so lutions comprise I conjugate 
pairs {(ii, j/i, jji, Zi)}. Recalling from above that 

the real and imaginary parts of each solution are linearly 
independent, no nontrivial solution may equal its com- 
plex conjugate (or even a scalar multiple thereof) so the 
solution pairs are all non-degenerate. Each solution pair 
yields a single multipole vector ±Ui via the cross product 
of its real and imaginary parts (Eqn. (|17|l '). 

Must this procedure find all the multipole vectors u^, 
as opposed to finding some several times each while 
missing others entirely? Yes. Pick an arbitrary but fixed 
Ui, and let (i, y, z) be a common solution of • (x, y, z) = 
and + + z"^ = 0. More concretely, any scalar 
multiple of {x, y, z) — {Uy+u1, —u^Uyzti u^, —u^Uz^i Uy) 
is such a solution, where {ux, Uy,Uz) arc the components 
of the given u^. A glance at H15|l shows that {x,y,z) 
must satisfy P(x,y, z) = as well, so our algorithm will 
find scalar multiples of {x,y,z) and (x,y,z) among the 
common solutions to P{x,y,z) = and x"^ + y'^ + z'^ — 
0. Substituting either {x, y, z) or (a;, y, z) into H17|l then 
gives the desired multipole vector u^, thus proving that 
our algorithm does indeed find all of them. 

In the non-generic case that some of the occur mul- 
tiple times each in (|15f) . does our algorithm get the mul- 
tiplicities right? The easy way to see that it does is to 
slightly perturb the polynomial P{x, y, z) so that the t 
multipole vectors become distinct. For the perturbed 
polynomial, exactly one pair of roots {(x, y, z), (i, y, z)} 
leads to each u^. If we gradually unperturb the polyno- 
mial, the \ii fuse back into their original groupings, and 
we see that each Uj gets mapped to by the right number 
of root pairs. 



III. IMPLEMENTATION 

The algorithm, implemented as a Math- 
ematica notebook, is freely available at 
|http: //www. geometrygames.org/Maxwell/ The user 
specifies a spherical harmonic as an array {a^o, ■ ■ • 
of coefficients in the usual spherical harmonic decom- 
position X]m=-£ '^fm^™- Note that the a^m are given 
explicitly only for m > 0; the values for m < the 
defined implicitly via ag^-m = (— 

As an example, the aim for the quadrupole component 
of Eriksen et al.'s Lagrange Internal Linear Combination 
(LILC) cleansing of the first-year WMAP data are 



,LILC 



{16.30,-2.57-1-4.98^,-18.80- 20.02i}. (18) 



Our Mathematica function AlmToPolynomial converts 
the 0^2^^ to a polynomial 



-19.67a;2-|-9.39y2-|-10.28z2 
-1-30.94 xy-f 7.71 yz -I- 3.98 zx. (19) 



P{x,y,z) 



The main algorithm, described in Section III Bl 
and implemented in the Mathematica function 





run time 


run time 


i 


without caching 


with caching 


8 


0.03 s 


0.02 s 


16 


0.2 s 


0.03 s 


32 


0.6 s 


0.1 s 


64 


3 s 


0.4 s 


128 


20 s 


2 s 



TABLE L The implementation of the PolynomialTo- 
MultipoleVectors algorithm runs quickly for low values of £. 
As (. increases, the call to the built-in Mathematica function 
Expand for simplifying a polynomial becomes a bottleneck 
(middle column). For Monte Carlo simulations, which re- 
quire many computations for the same £, pre-simplifying the 
spherical harmonics and caching them greatly improves per- 
formance (righthand column). 



PolynomialToMultipoleVectors, yields the Maxwell 
multipole vectors 



ui = { 0.975, 0.043, 0.220} 
U2 = { 0.634,-0.737,-0.234} 



(20) 



As an error check, one may feed the multipole 
vectors (|20|l to the function MultipoleVectorsTo- 
Polynomial to recover the harmonic polynomial H19|l . 



A. Speed 



The algorithm of Section FlI Bl as implemented in the 
function PolynomialToMultipoleVectors, runs exceed- 
ingly fast. On a 2.4 GHz personal computer, the algo- 
rithm computes the multipole vectors for randomly 
generated aem almost instantaneously for very small £, 
and within a fraction of a second for relevant larger £ 
(Tabled middle column). 

The algorithm's bottleneck lies in the built-in Math- 
ematica function Expand, with which PolynomialTo- 
MultipoleVectors reduces the polynomial in a (dis- 
cussed in Section III Bll to its canonical form cq + cia + 
C2a^ -I- • • • -I- C2ia^^ before solving for the roots. In 
particular, simplifying the polynomial to standard form 
is far slower than extracting the roots. When using 
Monte Carlo methods to compare the observed WMAP 
data to large numbers of simulated CMB skies, one may 
speed up the multipole computation by pre-computing 
the spherical harmonics, making the substitutions {x i— > 

— y 1^ 2a, zi— >i + 1 )} directly into the spherical 
harmonics, simplifying the polynomials to the canonical 
form bQ + bia + b2a^ + ■ ■ ■ + b2ecx^^, and caching the result. 
With pre-cached spherical harmonics, the runtime drops 
to a very manageable 0{£'^) (Tabled righthand column). 
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I 


2 4 8 16 32 64 128 


max Zuu' 


^Q_16 ;^Q-14 ^Q-ll 



TABLE II: The algorithm for computing the is robust. For 
each given value of I, the table shows the maximum angle (in 
radians) between a multipole vector Ui and the re-computed 



I 


23456 7 8 9 10 


maxZuDQTULiLc 


16° 3° 8° 6° 13° 16° 10° 12° 8° 



TABLE III: The method used to clean the WMAP data has 
a significant effect on the multipole vectors. For each value 
of I, the table gives the maximum angle (in degrees) between 
a multipole vector computed with the DQT data and the cor- 
responding multipole vector computed with the LILC data. 



B. Precision 

1. In theory 

The multipole algorithm is quite stable. To test it, 

1. begin with a random set of multipole vectors 
{ui, . . .,ue}, 

2. convert the multipole vectors {ui, . . . , u^} to a ho- 
mogeneous polynomial P, and then 

3. use P to recompute the multipole vectors 

{u'l, •■■,";,}. 

Table ^1 shovifs that the re-computed multipole vectors 
{u'l, . . . ,u'^} agree with the original multipole vectors 
{ui, . . . , u^} to high precision, even for fairly large values 
of ^. 



2. In practice 

Given the inherent numerical stability of the multi- 
pole algorithm, the real question becomes, how stable 
are the multipole vectors in a cosmological context? To 
answer that question, compute the low-i CMB multi- 
pole vectors for the first-year WMAP data, first using 
the DQ-corrected Tegmark (DQT) cleaning 0] and then 
using the Lagrange Internal Linear Combination (LILC) 
cleaning (J] of the same raw data. The resulting mul- 
tipole vectors differ by about 10° from one cleaning to 
the other (Table IIII(I . Thus we conclude that the com- 
puted multipole vectors have physical meaning, but only 
roughly. 

The most striking entry in Table ITTll is the 16° maxi- 
mum difference between the DQT and LILC quadrupole 
vectors. It's tempting to assume that this difference re- 
sults from the Doppler correction included in the DQT 
data but omitted from the LILC. However, the maximum 
difference between the non-Doppler-corrected Tegmark 
quadrupole vectors and the LILC quadrupole vectors is 



a whopping 37°. We must therefore concede that we 
know the true quadrupole vectors only very coarsely. 



IV. APPLICATIONS 

A. Quadrupole-octopole alignment 

Schwarz, Starkman, Huterer and Copi offer two vari- 
ations on their approach to measuring the quadrupole- 
octopole alignment 6]. Both variations begin with the 
two multipole vectors {u2,i, U2.2} for the quadrupole and 
the three multipole vectors {u3_i, U3_2, us ^} for the oc- 
topole. 

According to the unnormalized version, a cross product 
W2 = U2.1 X U2.2 defines a normal vector to the quadrupole 
plane. Similarly 

W3.I = U3^2 X U3,3 
W3,2 = U3^3 X U3,i 

W3,3 = U34 X U3,2 (21) 

define normal vectors to the three octopole planes. The 
three dot products Ai = |w2 • W3_i|, for i — 1,2,3, then 
measure the alignment of the quadrupole plane with each 
of the three octopole planes. Katz and Weeks take the 
sum S = Ai + A2 + A3, which they find to be unusually 
high at the 99.9% level when compared to 10^ Monte 
Carlo simulations . 

The normalized version is identical to the unnormal- 
ized one, except that it replaces the previous cross prod- 
ucts with unit length normal vectors = |w2| = 
|u2,i X U2,2| and w^j ^ = |w3,j|. 

So which version should we use? If our goal is to 
measure how well the quadrupole plane aligns with each 
octopole plane, then the normalized version works per- 
fectly. The unnormalized version, on the other hand, 
resists easy interpretation. By retaining the length of 
each cross product Uj j x Uj^v, it puts more weight on 
"well defined planes" (for which Uij and Uij/ are nearly 
orthogonal and which therefore robustly determine a nor- 
mal direction) and less weight on "poorly defined planes" 
(for which j and Uiji are nearly parallel or anti-parallel 
and which therefore only weakly determine a normal di- 
rection). Nevertheless it's not clear what overall signifi- 
cance the resulting statistic carries, nor whether this ver- 
sion is a good choice for the quadrupole-octopole compar- 
ison. For the moment let us remain neutral and compute 
statistics relative to both versions. 

The multipole algorithm of Section IIIBI as imple- 
mented in the Mathematica code of Section lTTll generated 
a million Monte Carlo simulations, whose quadrupole- 
octopole alignments were then compared to (a) the 
DQ-corrected Tegmark (DQT) cleaning of the first-year 
WMAP data and (b) the La grange Internal Linear 
Combination (LILC) cleaning of the first-year WMAP 
data jj]. Table IIVI shows the results. 
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unnormalized 
version 


normalized 
version 


DQT 


99.9% 


98.7% 


LILC 


97.3% 


98.7% 



TABLE IV: How unusual is the quadrupole-octopole align- 
ment observed in the first-year WMAP data? For each data 
set (DQT, LILC) and each version of the algorithm (nor- 
malizing the lengths of normal vectors or not), the table re- 
ports what fraction of a million Monte Carlo skies had weaker 
quadrupole-octopole correlations. 





unnormalized 
version 


normalized 
version 


DQT 


34.6% 


19.5% 


LILC 


17.9% 


21.0% 



TABLE V: Comparing the quadrupole {£ — 2) to the hexade- 
capole {£ = 4) finds the normalized statistic to be more stable 
than the unnormalized statistic, just as in the quadrupole- 
octopole comparison of Table HVI 



The normalized alignment statistic (Table IIVI right- 
hand column) appears far more robust than the unnor- 
malized alignment statistic (Table IIVI middle column) 
relative to small changes in the data, namely switching 
from the DQT cleaning of the first-year WMAP data to 
the LILC cleaning of the same data. Conceivably dumb 
luck alone might have given the DQT and LILC data 
similar values of the normalized alignment statistic. To 
test that possibility, one may compute the same statis- 
tics for other data. Comparing the quadrupole {£ — 2) 
to the hexadecapole {£ = 4) for the same data sets (Ta- 
ble again finds the normalized statistic more stable 
than the unnormalized one. Returning to our main data 
(Table I1V|) this suggests that the normalized statistic 
measures the quadrupole-octopole alignment more sta- 
bly than the unnormalized statistic does. 

The reason the unnormalized statistic rates the DQT 
data so much higher than the LILC data (Table HVI mid- 
dle column) has little to do with quadrupole-octopole 
alignment, but instead reflects the orthogonality of the 
two quadrupole vectors {u2,i,U2,2}- Specifically, their 
cross product W2 = U2,i x U2.2 has length |w2| = 0.990 
for the DQT data but only |w2| = 0.845 for the LILC 
data. So even if the LILC octopole bears exactly the 
same relation to the LILC quadrupole plane that the 
DQT octopole bears to the DQT quadrupole plane, the 
sum S = Ai + A2 + A3 will come out 17% higher for the 
DQT data than for the LILC data, because 



Wo 



0.990 
0.845 



1.17. 



(22) 



Indeed, for the actual data the ratio of the measured 
sums is 



DQT 



2.395 
2.048 



1.17. 



(23) 



Thus the near orthogonality of the quadrupole vectors 
{u2.i,U2.2} in the DQT data fully explains the differ- 
ence between the two entries in the middle column of Ta- 
ble ^3 The relative quality of the quadrupole-octopole 
alignments plays no role there. 

One concludes that the 99.9% value in Table |EI 
measures the combined unlikelihood of the quadrupole- 
octopole alignment and the near orthogonality of the 
DQT quadrupole vectors. To measure the quadrupole- 
octopole alignment alone, one should use the normal- 
ized statistic (Table llVl righthand column), which for the 
first-year WMAP data finds the alignment to be unusual 
at the 98.7% level. While this contradicts one recent 
claim 0, Q , it agrees well with an earlier estimate (using 
other methods) that the quadrupole-octopole alignment 
is unusual at the l-in-60 level 



B. Ecliptic plane alignment 

Schwarz et al. point out the alignment of the 
quadrupole and octopole vectors with the ecliptic plane. 
Here we test the statistical significance of that alignment 
and obtain results inconsistent with the 99.8% confidence 
level claimed in the preprint version of Ref. ^ but con- 
sistent with the revised ~ 99% confidence level appearing 
in the final published version of Ref. Q if one considers 
only a one-tailed distribution. 

First consider the normal to the quadrupole plane, 
which in galactic coordinates is (—106°, 57°) for the DQT 
data or (-112°, 62°) for the LILC data. Taking the dot 
product with the ecliptic pole ±(96.4°, 29.8°) gives 0.027 
or 0.080, nominally implying a confidence level of 97.3% 
or 92.0%, depending on the data set (Endnote 18]). 
However, assuming we would have been equally pleased 
with a dot product near 1 (meaning the quadrupole plane 
aligned with the galactic plane) , a two-tailed distribution 
is appropriate, dropping the confidence level to 94.6% or 
84.0%. 

Applying the same method to the three octopole planes 
gives dot products of 0.523, 0.045 and 0.179 (for the DQT 
data, in agreement with Ref. [^) or 0.555, 0.030 and 
0.146 (for the LILC data). At this point the Preprint Q 
computes a raw score which cannot be interpreted as a 
confidence level; to avoid that trap we take a different 
approach and examine the sum of the three dot prod- 
ucts, namely 0.747 (DQT) or 0.731 (LILC). Monte Carlo 
simulation of 10^ Gaussian random skies finds the sum 
to be larger than that 95.7% of the time (DQT) or 96.1% 
of the time (LILC). 

Finally, for more direct comparison with the results of 
Ref. '6'|, consider the joint sum of the single quadrupole 
dot product along with the three octopole dot prod- 
ucts. Monte Carlo simulation of 10^ Gaussian random 
skies finds this combined sum to be larger than the ob- 
served sum 99.0% of the time (DQT) or 98.9% of the 
time (LILC). This 99% result is weaker than the erro- 
neous 99.8% claim of Preprint but consistent with 
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4 


5 


6 


7 


8 


9 


10 


DQT 
LILC 


28% 
30% 


80% 
58% 


89% 
97% 


24% 
33% 


25% 
12% 


58% 
72% 


82% 
48% 


69% 
34% 



TABLE VI: For each data set (DQT, LILC) and each value of 
i (3. .10), the table shows what percentage of lO'' Monte Carlo 
simulations yielded a more planar set of multipole vectors. 



the 99.1% value in the final published Ref. |6|. Allow- 
ing for a two-tailed distribution — because quadrupole 
and octopole planes parallel to the ecliptic would have 
interested us just as much as ones perpendicular to the 
ecliptic — drops the confidence level to 98%, but even so 
the result remains valid at better than the 2ct level. 



C. Planarity of modes 

De Oliveira-Costa et al. found the WMAP octopole 
to be unusually planar at the l-in-20 level using their 
i-statistic |0. Multipole vectors offer an independent 
measure of planarity: the triple product \ det(ui, U2, U3)| 
of the three octopole vectors measures the volume of the 
parallelepiped they span. When the octopole vectors are 
mutually orthogonal, their triple product is one. At the 
other extreme, when the octopole vectors are coplanar 
their triple product is zero. Unfortunately the first-year 
WMAP octopole's triple product is only slightly low: 
28% (DQT) or 30% (LILC) of Gaussian random skies 
achieve a lower value. Thus we conclude that De Oliveira- 
Costa et al.'s t-statistic is more sensitive to the observed 
anomaly. 

Extending this approach to higher £, by summing the 



triple product of all three-element subsets of the multi- 
pole vectors {ui, . . . ,u^}, yields similarly disappointing 
results (Table lyijl . At £ — 5 the triple product method 
provides a weak 2(7 confirmation of Eriksen et. al's 3ct 
observation of non- planarity but otherwise the data 
in Table IVII seem consistent with a flat distribution on 
the interval [0%, 100%]. Again we conclude that the t- 
statistic better detects planarity or lack thereof. 



V. CONCLUSIONS 

Maxwell's multipole vector construction sheds much 
light on spherical harmonics and offers a clean and sim- 
ple approach to quantifying alignments among the low-£ 
CMB modes. In the case of the quadrupole-octopole 
alignment, simulations convincingly show the alignment 
to be unusual at the 98.7% level, corroborating earlier 
estimates of a l-in-60 alignment |0 while showing 
recent reports of 1-in-lOOO Q to have depended on 
the near orthogonality of the quadrupole vectors in 
the DQT data. The alignment of the quadrupole and 
octopole vectors with the ecliptic plane is confirmed at 
better than the 2<t level. 
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